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We study first-order spatial coherence for interacting polar bosons trapped in triple-well potentials. 
It is argued that besides the well-known coherence-producing couplings related to tunneling between 
the sites, there exists coherence based solely on intersite interactions, which prevails between the 
outer sites of the triple well when their total filling is odd. We find that it originates from the 
superposition of two degenerate many-body states in symmetric triple wells, and demonstrate its 
robustness against perturbations due to various tunneling mechanisms and pair exchange. 
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I. INTRODUCTION 

The study of the many-body physics of interact- 
ing particles stored in multiple-well geometries or peri- 
odic lattices has wide-ranging implications for the cross- 
fertilization of condensed-matter physics and quantum 
optics [![. Preparing ultracold polar molecules in such 
external scalar potentials possessing local minima, offers 
manifold opportunities to investigate the intricate phase 
diagrams related to the complex coherence properties im- 
plied by the long-range dipolar interactions between the 
molecules. In the simplest version of the Hubbard model 
with contact interactions only, the existence of first-order 
coherence is conventionally associated to the competition 
between kinetic (tunneling) energy and two-body interac- 
tions. By contrast, the existence of several energy scales 
associated with the nonlocal two-body interactions vastly 
enriches the structure of the phase diagram Q , and en- 
ables the creation of many phases with distinct coherence 
properties. 

For the observation of the effects of nonlocal interac- 
tions, a triple-well potential is among the conceptually 
simplest setups to present features which go beyond those 
of contact interactions. A major advantage of the triple 
well is that it is tractable numerically, e.g. for the ground 
state by exact diagonalization (ED), for a comparatively 
large number of bosons when the field operator expan- 
sion is restricted to the lowest-lying orbital per site, so 
that only three field-operator modes are involved. In op- 
tical lattice terminology, this would be referred to as the 
single-band Bose- Hubbard model [3[|. Previous theoreti- 
cal work has explored bosonic triple- well systems for con- 
tact interactions by a semiclassical approach mean- 
field treatments [H,Q, and the multiconfigurational time- 
dependent Hartree method for bosons [7|. Transistor- 
like effects in a triple well, where the occupation of the 
middle well controls tunneling between the outer wells 
were demonstrated in Q . Polar bosons in triple- well po- 
tentials, including however only the conventional single- 
particle tunneling mechanism, have been investigated by 
many-body Q and mean-field methods [Toj . 

In what follows, we demonstrate that a triple well stor- 
ing polar bosons with dipolar interactions can exhibit 



variants of the first-order coherence of matter waves of 
physically distinct origins. Besides the coherence due to 
single-particle and occupation-number-weighed tunnel- 
ing we show that, in the appropriate parameter regime, 
the intersite interaction can generate macroscopic first- 
order coherence between the two outer sites of the triple 
well, in the absence of any kind of tunneling. We demon- 
strate both the origin of this interaction-induced co- 
herence and its stability against perturbations due to 
various tunneling mechanisms (both conventional short- 
range and long-range) as well as pair exchange. We also 
point out that a possible destabilization by a possible 
asymmetry of the trap potential can be counteracted by a 
small amount of occupation-number-weighed tunneling, 
which stems from interactions, and analyze the finite- 
temperature behavior of the interaction-induced coher- 
ence and propose its possible measurement. Finally, we 
discuss specific features of the phase diagram of the triple 
well related to pair exchange and occupation-number- 
weighed tunneling, establishing inter alia a procedure to 
measure the latter. 



II. TRIPLE WELL HAMILTONIAN 



We consider N dipolar bosons trapped in triple-well 
potentials where both the potential depth in each well 
and the relative orientation of the polarized dipoles is 
adjustable [|, [lo[. It is furthermore assumed that the 
lowest-energy orbital wavefunction in each well, 4>i(Y), 
where i denotes the well index, is fixed and independent 
of the particle numbers per site (the filling). The dipole 
moment of the bosons is sufficiently strong in the sense 
that a large dipole coupling dominates the contact inter- 
action coupling; this is realistically obtained with ultra- 
cold polar molecules with large permanent electric dipole 
moment [Til ] . Expressing the bosonic field operator as 
^( r ) = £i=i 4>i( r )0'ii the low-energy Hamiltonian for the 
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triple well is of the single-band Bosc-Hubbard form 
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Here ^ = -/dr#(r) [-V 2 /2m + Vt rap (r)] & +1 (r) 
(with Vt ra p the triple well scalar potential and we set 
H = 1), J 2 = -f\Mr)\ 2 tf(r')&+i(r')V dd (r - r')drdr', 
and W = /^(r)^(r')^ + i(r)^ +1 (r')Kid(r-r')drdr' 
are the single-particle, occupation-number-weighed, and 
pair-tunneling rates respectively [HI, EH- The cou- 
plings U = gf\Mr)\ 4 dr + f |^(r)| 2 |^(r')| 2 Vdd(r - 
r')drdr' characterizes the on-site interaction, and V = 
J |0i(r)| 2 |(/> i+ i(r')| 2 Vdd(r - r')drrfr' is due to dipolar- 
interaction-induced coupling between sites; we neglect 
the contribution of the contact contribution to the 
J2,V,W terms. Finally, e represents identical single- 
particle energies at each site and 8E denotes a bias energy 
applied to site 1, while g = 47ra s /ra is the short-range in- 
teraction constant with a s the s-wave scattering length. 
The term X)»=i e "-i = ^ s displayed for completeness 
here, even though it adds only a constant to the Hamilto- 
nian: We choose e as the (arbitrary) energy scale below. 

The dipolar interaction between polarized dipoles at 
positions r, r' is given by Vdd(r— r') = d 2 (l — 3cos 2 6>)/|r— 
r'| 3 , where 9 is the angle between r — r', and the dipolc 
moment d. d 2 = /j,Qfi 2 / (Air) and d 2 /(Aireo) for magnetic 
and electric dipoles, respectively, where fi is the magnetic 
dipole moment and d is the electric dipole moment. From 
point-like, tightly localized orbitals <fii(r) to two infinite 
chains of dipoles where the width of two aligned orbitals 
is much larger than its separation (e.g., in Fig. [S] below, 
the condition o~ y 3> er z amounts to two infinite chains 
of dipoles), the factor a for intersite interactions varies 
from 8 to 4 [{J . The ratios a of next-nearest-neighbor and 
nearest-neighbor sites vary in general somewhat for the 
contributions J2, W, and V, respectively, depending on 
the geometry of the triple- well potentials and the overlap 
of orbital wavefunctions. For simplicity and because it 
does not affect our main conclusions below, we take their 
ratio factors a to be equal; we choose a = 8. 

We include the pair tunneling W, which has been 
shown to decide on the question of coherent versus frag- 
mented many-body states in general two-mode models 
[l4l ]. On an optical lattice, it was subsequently demon- 
strated that pair-tunneling can change significantly the 
phase diagram of polar bosons in the extended Bosc- 
Hubbard model [13[. 

The fundamental object of interest in what follows, the 
many-body wavefunction, can generally be expanded in 



a Fock basis, 

I*) = X] /( n i! n 2,n3)|ni,n2,n 3 ), (2) 

111,112,113 

where n\ + rii + n 3 = N. 

III. VARIANTS OF COHERENCE 

A. Vanishing tunneling couplings 

In the absence of any tunneling, i.e., Ji = J2 = W = 0, 
the total energy of the triple-well system subject to the 
Hamiltonian ([!]) reads, 



£(ni,n 3 ) = (V - U)(N - m - n 3 )[m + n 3 ] 
+ (V/a - U)n in3 . 



(3) 



We concisely review the phase diagram of polar bosons 
in a triple well as found in Phase A, cf. Fig.[TJ is 
characterized by m = n 3 , n 2 ^ 0, and appears for U > 
and V < j^U, and U < and V < SU. Phase B occurs 
for m = n 3 = f , n 2 = 0, for f^U < V < SU and 
U > 0. The case that V > U for U < 0, as well as 
V > 8U for U > leads to the phase C where n\ = 
n 3 = 0. Finally, for phase D macroscopic states with 
111 = and n 3 ^ and n\ 7^ and n 3 = are coherently 
superposed in a Schrodinger cat type state, where the 
regime of parameters is U < and 8U < V < U. 
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FIG. 1. (color online) (a) Phase diagram in terms of the 
magnitude of n/N = (hi + h^,)/N as a function of U and V 
for iV = 30, J\ = J2 = W = 0. The dashed lines in panel (a) 
show the boundaries between the phases A-D that are shown 
schematically in (b). 

We now point out that in Ref. , an important prop- 
erty of phase A has been overlooked. For phase A, in 
which n 2 is restricted to be integer by the fact that this 
minimizes the energy Eq. @, both ri\ and n 3 can be frac- 
tional (and not only integer, as stated in 0). We found, 
as also numerically verified by ED, that when the average 
fillings in each site, n% = (n,) fulfill n = ni + n 3 is odd, 
so that fi\ and n 3 are fractional there exist two-fold de- 
generate states in phase A, with the degenerate subspace 
spanned by the two states 




|*i) = |(n-l)/2,JV-Ti,(n + l)/2), 
|*2> = |(n+l)/2,JV-Ti,(n-l)/2). 



(4) 
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On the other hand, when n is even, the many-body state 
is a single Fock state, \V) = \n/2,N - n,n/2). All off- 
diagonal first-order coherence matrix elements are then 

zero, i.e., = 923' = 9 < i3 = 0- This is, however, not 
the case for a superposition of the states in Eq. Q , as we 
will now argue. 




FIG. 2. (color online) Distribution of many-body amplitudes 
in the wave function ([2]) for N — 30 particles and odd n — 17 
(left, U = 2.5, V = -3, N = 30, n = 17) and even n = 18 
(right, U = 4, V = —3). The switching between the double- 
and single-peak structures corresponds to the "coherence fan" 
structure in shown where macroscopic coherence between 
sites 1 and 3 is obtained only for odd n. The double-peak 
structure on the left is robustly obtained also when small 
tunneling couplings are turned on, cf. Fig.[3](b)-(d). 

The many-body state is in general any linear combina- 
tion of two degenerate states, i.e., \^f) — a\^>i) + (3\^2) 
where the coefficients a and /? can be any complex num- 
ber and fulfill \a\ 2 + \(3\ 2 = 1. However, for any small 
tunneling in a realistic system, e.g., any small J\, J 2 , or 
W, a and /3 adjust to relative phase zero. Moreover, the 
spatial symmetry between sites 1 and 3 requires fi\ = 7x3. 
The only possible form of the many-body wavefunction 
then is |\fr) = ^=[|*i) + *2>], with some global phase 0. 
Setting the latter to zero, we conclude that 



l*) = -^Pi> 



1*2 



(5) 



The above many-body state induces macroscopically 
large nonlocal first-order coherence, i.e., the quantity 
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is of order unity, cf. Fig.[3j where n = 5Zi=i 

We stress that the interaction-induced coherence only 
exists between the outer sites of the triple well, i.e., 

012 = 923 = and 9i3 = TT' The single-particle 
density matrix indeed has three macroscopic eigenvalues, 
which are f , ^j^ 1 , and N—n. This demonstrates that, 
seen from the point of view of the entire triple well, the 
superposed states with macroscopic interaction-induced 
coherence between the outer sites, are globally speaking 
fragmented condensate states [l5j . 

We point out that the superposition (|5|) of two-fold 
degenerate states is much more robust than the degener- 
ate states occurring in phase C. For vanishing tunneling, 
the many-body wavefunction in the parameter regime of 
phase C should be any of three degenerate states, i.e., 



I*) = a\N, 0, 0) + 6|0, N, 0) + c\0, 0, N) where any two of 
a, 6, and c are zero and the remaining one unity. But this 
tri-fold degeneracy is broken readily against a small per- 
turbation of first-order coupling, e.g., the single-particle 
tunneling between wells, and the many-body wavefunc- 
tion favors the single-particle state |0, N, 0). By contrast, 
a small perturbation due to tunneling does not destroy 
the degeneracy of |\E'i) and ^2), as we will demonstrate 
in the following subsection. 



B. Finite tunneling couplings 
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FIG. 3. (color online) Illustration of the insensitivity of the 
"fan" structure of the coherence function g[]^ , related to the 
phase transitions between even and odd values of n for varying 
tunneling couplings. Magnitude of 3^3' f° r ( a ) no tunneling, 
Ji = J 2 = W = 0, (b) Ji = 0.02, and J 2 = W = 0, (c) J 2 = 
0.001, and Ji = W = 0, (d) W = -0.02, and Ji = J 2 = 0; 
N - 30 in all plots. 

The emergence of the first-order coherence for n = odd 
and disappearance for n = even results in the formation 
of a fan structure of the first-order coherence in the pa- 
rameter regime of phase A, see Fig(3] To explain the 
structure of the phase diagram in region A of the phase 
diagram in more detail, we now argue that while a small 
amount of single-particle, occupation-number-weighed, 
or pair tunneling slightly broadens the statistical distri- 
bution of the many-body wavefunction ([2]) around the 
superposed state in Fock space, it does not change the 
structure of the phase diagram. For concreteness, with a 
finite but small Ji, the many-body wavefunction should 
be of the form, 

\*) = Cl \(n-l)/2,N-n,(n+l)/2) 
+ c 2 \(n+l)/2,N-n,(n-l)/2) 
+ c 3 |(n-3)/2,iV-n,(n + 3)/2) (7) 
+ c 4 |(n + 3)/2, N-n,(n- 3)/2) 
+ other terms . . . , 
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where the additional single-particle states are determined 
by Ji and the larger J 1; the more single-particle states 
need to be included. For small J\, c\ and C2 are 
dominant, and determined by U and V, such that the 
interaction-induced coherence is obtained. Only with 
strongly increasing Ji, the contributions of more single- 
particle states broaden the distribution of the many-body 
state in Fock space, significantly reducing the value of c\ 
and C2 , and then suppressing the interaction-induced co- 
herence. A similar discussion can be applied to the cou- 
plings J 2 and W. Therefore, effects caused by the various 
tunnelings do not result in an essential change for the 
interaction-induced coherence, as long as the tunneling 
perturbations remain sufficiently small. This robustness 
of the coherence against perturbations has been numeri- 
cally verified by the the ED calculations shown in Fig. [3] 
It can be concluded from Fig. [3] that the magnitude of 
the coherence changes due to finite tunnelings of various 
origin. However, the phase boundaries (the "fan" struc- 
ture), while slightly blurred by small tunneling, remain 
basically intact. Note for the proper interpretation of the 
parameter values chosen that the relevant quantities to 
compare with each other are J\ and nJi- 

By comparing Figs.|3](b) and (c) we see that the nature 
of the crossover between the phases B and C is strongly 
influenced by number- weighed tunneling, while ordinary 
tunneling has basically no effect on the transition. This 
highlights the long-range nature of the tunneling mech- 
anism correlated to site occupation, embodied in the J2 
terms. In addition, it furnishes a possible method to 
measure J2 independently from J\ in the triple well by 
investigating the transition region between phases B and 
C. 



We finally conclude from Fig.[3](d) that the effect of 
pair-exchanges oc W on the single-particle coherence 
magnitude is most pronounced, confirming the profound 
effect pair exchanges have on the first-order coherence 
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C. Small bias energies 

When a small bias energy 8E (chosen here to be nega- 
tive) is accounted for, the actual triple-well ground state 
should be | ^J^, N — n, and the associated energy 

is (V - U)(N - n)n + (£-{/) "ipl + SE(n + ±) while 
the energy of the superposition state with respect to the 
Hamiltonian ©) is (V -U)(N- n)n + - U) r -^f± + 
6En (the superposition (|5|) is not the eigenstate of the 
Hamiltonian ([!])). Therefore the superposition energy 
is higher than the ground-state energy by \SE\/2, indi- 
cating that for any asymmetry between sites 1 and 3, 
the two-fold degeneracy and therefore the coherence is 
broken. However, for a negative intcrsite interaction, 
i.e., V < 0, choosing Gaussian orbitals (also see sec- 
tion IIVI below) shows that J2 is positive. This implies 
that J2 prefers the superposed states and stabilizes the 
interaction-induced coherence states. For a small 5E < 
and J2 > 0, the ansatz for the many-body state can be 
written 77 1 2±I, N - n, ^)+7 \^,N-n, s±±), where 
77 and 7 are real and have equal sign, with normalization 
rf + 7 2 = 1. The energy contributed by the bias SE is 
[f + i(?? 2 - 7 2 )] 5E, while that of J 2 is -^77777(77, + 1). 
It is then straightforward to obtain that the condition 
for stabilizing the interaction-induced coherence state is 
given by J 2 > for not too large |<5-E| (for increasing 

J 2 , more single-particle states need to be included). This 
condition has been numerically confirmed as well. 

D. Detection of interaction-induced coherence 

We propose to verify the existence of purely 
interaction-based coherence with the nonequilibrium dy- 
namics of the nonlocal coherent state. The initially as- 
sumed superposition state including the ground-state co- 
herence is quenched towards a nonequilibrium state by 
breaking suddenly the symmetry of sites 1 and 3, by ap- 
plying a single-energy site bias to site 1 instantaneously, 
cf. the illustration in Fig.|4j The time-dependent many- 
body wavefunction is then of the form, 



l*(*)) = 77= exp 



it{ (V - U){Nn - n z ) + 



a J 2 



Ne 



^cxp 
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(8) 



where n is assumed to be odd. Therefore, it is obtained 
readily that 

= ^H±2l cos [SEt] (9) 
in 

9^=9^=0- (10) 



Formulas ([9]) and pop show that there exists exclusively 
the time-dependent phase coherence between sites I and 
3, for which the out-of-equilibrium oscillation frequency 
is determined by 5E, in complete agreement with our 
numerical results displayed in Figj4] The interaction- 
induced coherence therefore can be measured by means 
of atomic interference, where at different instants of evo- 
lution with respect to the same initial setup, the trap po- 
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tcntials is released such that the freely expanding clouds 
from sites 1 and 3 will interfere. The visibility of the 
interference pattern then amounts to a measure of the 
first-order coherence [l6j . 
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FIG. 4. (color online) Time evolution of (t) in the absence 
of tunneling, Ji = J2 = W — 0, for an initial state with 
U = 1.5, V = —2, and 5_E = 0, quenched into 8E = —1 (see 
schematic representation on top) .The black solid line is from 
a numerical calculation, and the red dashed line from formula 
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E. Coherence versus temperature 

Here, we explore the effect of temperature on the 
interaction-induced coherence. In the canonical ensem- 
ble, the thermal average of an operator O is (O) = 

N 

J2j=o e z i^j \O\Ej), where the canonical partition sum 



Z = Yjj=o e Ej l T - In the zero temperature limit, (O) 



yields the expectation value of O in the ground state. 
Fig [5] shows that in the regime of small temperatures, 
e.g., T < 0.5, g[]^(T) is only slightly reduced, indicating 
that the low-energy excitation states do leave the super- 
posed states basically intact. With increasing T > 0.5, 
the high-level excitation states yield to a significant re- 
duction of g[^(T). 
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FIG. 5. (color online) The first-order spatial coherence, 
9ij (T) as a function of temperature, with U = 2.5, V = —3, 
N = 30, and Ji = J 2 = W = 0. 



We suggest to realize the interaction-induced coherence 
in small samples because large atom numbers suppresses 
the fan structure window size shown in Fig. [3] 

A triple-pancake trap 2D configuration to achieve U > 
and V < is shown in FiglHl To understand the pa- 
rameter relation between V, J2, and W, we take as the 
low-energy approximations to the exact Wannier states 
orbital wavefunctions correspondin g t o the harmonic- 
oscillator ground-states at each well [18| . Specifically the 
single-particle orbitals arc assumed to be of the form, 



Mr) = 
fair) = 



1 
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1 



= exp 



exp 



exp 



Aal 4<r2 
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4a 2 z 
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4cr 2 



6(x), (11) 



S(x). 



We take the separation of the two nearest-neighbor or- 
bital wavefunctions, I, as the unit of length and a z to ob- 
tain the relation of —J2/V and W/V as shown in Fig. [7] 
(the ratios are independent of a y ). The numerical results 
upon evaluating the parameter ratios show that J2 > 
and W < are consistently smaller than \V\, and with 



IV. EXPERIMENTAL SETUP AND 
PARAMETER CORRELATIONS 

We now discuss experimental feasibility and realistic 
parameter values for dipolar interactions. Ref. pro- 
vided a rather detailed scheme of how to realize the re- 
quired triple well configuration experimentally In partic- 
ular, using an acousto-optic modulator to toggle the dim- 
ple formed by optical potentials between several positions 
at high rate to create almost arbitrary time-averaged po- 
tentials [TtJ can be applied to our case to adjust the 
single-energy offsets. It should not be difficult for exist- 
ing experiments to reach the parameter regime of phase 
A which is quite broad in terms of the values of U and V. 
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FIG. 6. (color online) The configuration of a triple pancake 
system where the three 2D pancakes are aligned along the y 
direction, and all dipoles are oriented along z (red arrows), 
such that U > and V < 0. 
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FIG. 7. (color online) The ratios —J2/V and W/V for vari- 
ous values of the Gaussian width <r z , cf. Fig[6] Filled squares 
(—J2/V) and circles (W/V) are numerically evaluated data 
points, and lines are a guide to the eye. 

increasing overlap of two orbital wavef unctions, the ra- 
tios —J2/V and W/V increase monotonically (see Fig.[7|). 
By contrast, for contact interaction, W/V is constant and 
— J2/V 3> 1 for small overlap, decreasing with increasing 
overlap [l||. This indicates that the dipole-dipole inter- 
action tends to enhance the interaction-induced number- 
weighed tunneling over the nearest-neighbor interaction 
coupling for increasing overlap (in the setup of Fig.[7J) , 
while contact interactions have the opposite effect. 

V. CONCLUSION 

We have shown that intersite interaction in a triple 
well, due to large dipolc-dipolc interactions, in a broad 
parameter regime, can generate macroscopic first-order 
coherence between the two outer sites, in the absence 
of tunneling between any of the three sites. This 
interaction-induced coherence originates from the super- 
position of two degenerate states of polar bosons in 
a symmetric triple-well potential. While globally the 
many-body state is fragmented, it exhibits this par- 
tial coherence between the modes localized on the outer 



sites, which is robust against tunneling perturbations and 
against pair exchanges. We have also demonstrated that 
the fragility of the coherence against site biases (a mis- 
match of energy minima lifting the degeneracy) can be 
counterbalanced by a small contribution of occupation- 
number-weighed tunneling. 

The interaction-induced coherence between sites 1 and 
3 executes Josephson oscillations when the initial super- 
position state is quenched towards a nonequilibrium state 
by suddenly breaking the symmetry of the triple well 
(rapidly turning on a site bias) , where the associated os- 
cillation frequency is determined by the bias. This pro- 
vides an experimental method to measure the strength 
of the the interaction-induced mesoscopic coherence be- 
tween the outer sites of the triple well. 

Finally, we note that there exists a recent work in 
which it was shown that repulsive interactions induce 
first-order correlation between spatially separated chains 
of one-dimensional dipolar fcrmions [19( . Our origin of 
coherence is determined by both the spatial symmetry of 
the triple well and the intersite interaction. On the other 
hand, due to the existence of many degenerate states in 
a very narrow window of intertube interactions in the 
fermionic system, an infinitely small intertube tunneling 
can mix many degenerate states, and then lead to a co- 
herent ground state. Although the nonlocal interaction 
coupling V is the crucial factor in both cases, the in- 
teraction in the fermionic case is repulsive, i.e., V > 0, 
in distinction to the present case of elementary bosons, 
where the first-order coherence arises from negative V. 
In addition, the coherent-state regime for V is very nar- 
row in the fermionic system, in contrast to the bosonic 
triple well, where V and U cover a wide area in parameter 
space. 
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